function a=DLS(x,y)
G = normalequation(x);
c = RHS(x,y);
a = G \ c;